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The three-band model relevant to high temperature copper-oxide superconductors is solved us- 
ing single-site dynamical mean field theory and a tight-binding parametrization of the copper and 
oxygen bands. For a band filling of one hole per unit cell the metal/charge-transfer-insulator phase 
diagram is determined. The electron spectral function, optical conductivity and quasiparticle mass 
enhancement are computed as functions of electron and hole doping for parameters such that the 
corresponding to the paramagnetic metal and charge-transfer insulator sides of the one hole per cell 
phase diagram. The optical conductivity is computed using the Peierls phase approximation for the 
optical matrix elements. The calculation includes the physics of "Zhang- Rice singlets" . The effects 
of antiferromagnetism on the magnitude of the gap and the relation between correlation strength 
and doping-induced changes in state density are determined. Three band and one band models 
are compared. The two models are found to yield quantitatively consistent results for all ener- 
gies less than about 4eV, including energies in the vicinity of the charge-transfer gap. Parameters 
on the insulating side of the metal/charge-transfer insulator phase boundary lead to gaps which 
are too large and near-gap conductivities which are too small relative to data. The results place 
the cuprates clearly in the intermediate correlation regime, on the paramagnetic metal side of the 
metal/charge-transfer insulator phase boundary. 



PACS numbers: 71.27.+a, 71.10.Hf, 71.30.-|-h, 74.72.-h 

I. INTRODUCTION 

The interplay of band theoretic issues of hybridization 
and chemical bonding with the quantum chemical issue of 
strong local correlations is basic to the physics of many 
important materials. In an important paper^ Zaanen, 
Sawatzky and Allen classified insulating transition metal 
oxides as "charge-transfer" or "Mott" insulators accord- 
ing to whether the physics could be discussed solely in 
terms of strong correlated transition metal d states or 
whether transitions to O 2p states were important to the 
low energy physics. The issue arises with particular force 
in the high temperature superconductors where the Cu 
(f and d}^ states are not far in energy from O 2p states, 
but the Cu state is very far away in energy. In this cir- 
cumstance a strong particle- hole asymmetry is expected^ , 
with doped electrons residing on Cu sites whereas doped 
holes reside mainly on the O, but may be bound to Cu 
spins creating "Zhang-Rice singlets" ^ 

Quantifying this appealing physical picture requires 



solving an electronic structure problem with multiple 
scales, including a correlation energy on the Cu site 
~ 8 — lOeV, a Cu-0 energy level difference of 
and a Cu-0 hybridization ~ 1.6eVjf In this paper we 
use single-site dynamical mean-field theory (DMFT) to 
solve a model involving both copper and oxygen orbitals, 
developing a comprehensive theoretical picture of the 
electronic structure and optical conductivity of undoped 
and doped cuprate materials across the charge-transfer- 
insulator to charge-transfer metal phase diagram, includ- 
ing the effect of antiferromagnetism on the spectra and 
optics. We use newly improved Exact Diagonalization 
(ED)^'^ and Continuous-Time Quantum Monte Carlo 
(CT-QMC)S impurity solvers. These methods have dif- 
ferent sources of error and we find consistent results with 
the two methods. 

Our work is related to previous work of Dopf et al.^^ 
Georges et al.r^ and Zolfl et ali^ who each studied 
one particular parameter value. Modern developments 
in computers and solvers mean that we are able to ob- 
tain much more information. Our work also has some 
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overlap with more recent work of Macridin et aiJ^ who 
used the Dynamical Cluster Approximation on a four site 
cluster to study momentum dependence and the onset of 
superconductivity. We compare our findings to very re- 
cent work of Weber et ali^ who used similar methods to 
study a similar model. Weber et al. focused on specific 
parameters; we focus on the spectral functions and con- 
ductivities over a wider energy range, varying the charge- 
transfer gap to explore all regions of the theoretical phase 
diagram and present a comparison of the low energy be- 
havior of the copper-oxygen model to that of an effective 
one-band model. 

We present evidence that a one-band model provides a 
reliable picture of the spectral functions and conductivity 
for frequencies less than about 4eV (note that this range 
extends about a factor of two above the charge-transfer 
gap in frequency). We find an electron- hole asymme- 
try in the self-energy. The asymmetry is much more 
pronounced for parameters such that the undoped ma- 
terial is a charge-transfer insulator (insulating even in 
the paramagnetic phase). The asymmetry, however, is 
not reflected in the Fermi velocity renormalization or the 
low frequency optical matrix oscillator strength, where 
changes in the electronic structure compensate for the 
differences in correlation strength. Comparison of our re- 
sults to data suggests that the cuprates are on the metal- 
lic side of the single-site DMFT phase diagram, with an- 
tiferromagnetism being essential to produce the gap in 
the undoped material. 

The rest of this paper is organized as following: in sec- 
tion [IT] we describe the model, in section IIIII we present 
the calculated phase diagram, spectral functions, self- 
energies and optical conductivities, in section lTVl we com- 
pare the copper-oxygen model results to those obtained 
from computations performed on the one-band Hubbard 
model. Section IVl contains a summary of our results and 
a conclusion. 



II. MODEL 

We analyze the canonical two dimensional "copper- 
oxygen" Hamiltonian^'^'^ retaining the Cu d.j;2_y2 and 
O 2pa- orbitals (whose momentum (p) components are 
created by the operators djj^,plp^,p'^yp^). We allow 
for the possibility of two-sublattice antiferromagnetism 
by doubling the unit cell. The Hamiltonian is there- 
fore a 6-band model H = i/eband + Hmt in the mag- 
netic Brillouin zone. To write the band theoretical 
part we divide the lattice into two sublattices, A and 
B, distinguish the oxygen sites displaced from the Cu 
in the x and y directions, adopt the basis {ip) — 

i'^Apa'PA,x,pa'PA,y,pa^'^Bpa^PB,x,pa^PB,y,p<T) WritC 



where the 3x3 matrices are 
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and Hint — U J2i '^dt^di- Here we neglect oxygen-oxygen 
hopping. We use the value tpd = 1.6eV suggested by 
band theory calculations.^ 

Because only the Cu site is interacting, we may in- 
tegrate out the oxygen band to obtain an effective one- 
orbital model, which we solve in the single-site DMFT— 
using EDl>^ and CT-QMG^ methods described in the lit- 
erature. In the ED calculations typically 9 bath states 
were used and the results were verified by occasional 
large calculations; for CT-QMC temperatures studied 
were typically T — tpd /2b for the phase diagram and 
T = tpd/16 for spectral functions and conductivities. 
The main result of the DMFT calculation is a self-energy 
which, in the antiferromagnetic phase, is spin dependent 
and takes the form 
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with Sa,t — S-B,l- 

In the single-site DMFT method the self-energy is mo- 
mentum independent so the conductivity may be com- 
puted fromi^ 

X Tr [j6band(p)A(tJ + f2,p)j6band(p)A(a),p)] , 

where the integral is over the magnetic Brillouin zone, 
the spectral function A is 



A(t^,p) = -{G{uj,T.{uj + ie,p)-G(uj,T.(uj-ie,p)), 

(6) 

the matrix Green function G at frequency z and chemical 
potential /i is G{z,p) = (zl + ^ — S(z) — i?6band) ^, the 
current operator j — 5'ii/5px is 



J6band ~ ^ ^ 



(1) 






4 Px 

e 2 








e 


o\ 


e 2 








— e 2 





























e 2 








e 2 





-e 2 































0/ 



(7) 



3 



11 



10 



> 

CD 



QMC - 
















EC 


) ■■ 










/, J 




- 

■ 


















- 
- 


Insulator 










- 
- 










// ^ 






- 






























i/letal 




















- 


3 


-6 - 


4 - 


2 ( 


) 'c 




i 6 8 



A=£p-£jj [eV] 

FIG. 1: Metal-insulator phase diagram in plane of interaction 
U and p — d level splitting A for one hole per Cu02 unit 
in paramagnetic phase. Dotted lines (blue online): phase 
boundaries from ED calculation at T = ; solid lines (red 
online) indicate limit of stability of metallic phase from CT- 
QMC calculation at T = l/40eV. In the region between the 
two lines metallic and insulating solutions coexist. Square 
(blue online), circle (black online) and triangle (red online): 
parameters studied by Ref. [lol Ref. Ref. |l^ respectively. 



and f{uj) is the Fermi function. 

For comparison we have studied the one-band model 



described by H 
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t — 0.37eV is chosen to reproduce the non-interacting 
bandwidth of the Cu-0 antibonding band passing 
through the Fermi level and UcS fixed so as to reproduce 
the correlation gap. 



III. RESULTS 



A. Phase diagram 
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FIG. 2: Electron spectral function per spin for paramag- 
netic three-band model (negative frequency: removal spec- 
trum; positive frequency: addition spectrum) (solid lines) and 
projections onto d (dashed; red online) and p (dotted line; 
blue online) states calculated with the ED solver. Parame- 
ters: U = 9eV and A = 4eV (upper panels) and A = 2eV 
(lower panels). Upper graphs: 0.15 hole doping (A — 4eV: 
Ed = -7.7eV, £p = -3.7eV; A = 2eV, Sd = -6.3eV, Sp = 
— 4.3eV) middle graphs: undoped (A = 4eV, Ed ~ — 8.8eV 



-4.8eV and A = 2eV, Ed 



-7.9eV Er, 



-5.9eV). 



Lower graphs: 0.15 electron doping (A = 4eV, Ed = — 9.8eV, 
Ep = -5.8eV; A = 2eV; Ed = -9.6eV, Ep = -7.6eV). 



Fig. [T] shows the boundary between metallic and para- 
magnetic insulating solutions, calculated for one hole per 
Cu02 unit as a function of interaction strength U and 
Cu-0 energy level splitting A(= Sp — Sd) in the paramag- 
netic phase. As in the single-site DMFT of the one-band 
Hubbard modeli^ a coexistence region is observed where 
both metallic and insulating solutions exist. We define 
the insulating and metallic boundaries of the coexistence 
region at fixed large U to be Ac2 and Ad respectively. 
The phase diagram is obtained at T = using the ED 
solver; the results were verified using CT-QMC by scan- 
ning A at selected U values. We find almost perfect 
agreement for Aci; the CT-QMC calculation of Ac2 line 



represents the smallest A value at which a metallic phase 
can be found at temperature T = l/40eV. The known 
T-dependence of the single-site DMFT boundary for the 
single-band Hubbard model^^ suggests that the discrep- 
ancy may simply be a finite-temperature effect related to 
the different temperatures used in the two calculations. 

B. Spectral function 

Fig. [^Ipresents the many body density of states (DOS) 
(electron removal spectrum for energy < and electron 
addition spectrum for energy > 0) as well as the pro- 
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jections onto the Cu-d and 0-2^0- states calculated in 
the paramagnetic phase using the ED method. The zero 
of energy is the chemical potential. The upper panel 
presents results for A > Ac2 (so within single-site DMFT 
the paramagnetic phase of the undoped material is metal- 
lic) while the lower panel presents results for A ~ Ac2 
such that the undoped material is a charge-transfer insu- 
lator. The main difference between the two spectra is the 
presence or absence of a gap in the undoped material. 

At lowest energy (binding energy ^ lleV) a peak is 
seen, of mainly c?-character. This peak corresponds to 
removing one electron and leaving the Cu in the state; 
it is pushed down from the bare energy by a level 
repulsion due to hybridization with the O states. Thus it 
is not correct to identify the position of this peak directly 
with U, as is sometimes done in literature. 

In the binding energy range 4eV ^ s < 8eV a mainly 
oxygen-like band is seen. The very sharp peak corre- 
sponds to the non-bonding oxygen state; it would be 
broadened if oxygen-oxygen hopping were taken into ac- 
count (although the vanishing of the Cu-0 hybridiza- 
tion at the r point means that a singularity would re- 
main at the bare oxygen energy). The broad structure of 
mixed oxygen and copper character lying below the non- 
bonding state may be thought of as the "bonding" linear 
combination of Cu and O states pushed down below the 
non-bonding O level by hybridization to the Cu d level. 
This feature was identified in Ref. [ij as the copper up per 
Hubbard band corresponding to the d^ state, but Ref. [Ij 
did not present results over a wide enough range to de- 
termine if the ^ lleV peak (which we find to correspond 
to the d^ state) was present in their calculations. 

The structure of mixed Cu-0 character at binding en- 
ergies in the range e ^ 1 — 4eV corresponds to the Zhang- 
Rice singlet states. Calculations (not shown here) in 
which the Cu orbital is forced to be fully spin polar- 
ized show that these states correspond to holes with the 
same spin as the deep-lying removal state. The states 
in the electron addition spectrum are mainly of Cu c?^" 
character, and play the role of the upper Hubbard band. 

The effect of antiferromagnetism on the magnitude of 
the gap in the charge-transfer iiisulator state has been 
the subject of debate, with Ref. 18 arguing on the basis 
of Hubbard model calculations that antiferromagnetism 
increases the gap significantly while the conductivity cal- 
culations presented in Ref. 14 were interpreted as indicat- 
ing no significant effect of antiferromagnetism on the gap. 
We have used the "quasiparticle equation" method of 
Ref. to determine the gap values a,tU — 9eV, A = 2eV 
finding that the gap in the paramagnetic insulating phase 
is 2.87eV while the addition of antiferromagnetism shifts 
the gap to 3.47eV. Some of the difference between our re- 
sults and those of Ref. H may arise from the extremely 
small value of the calculated paramagnetic-state conduc- 
tivity in the near gap region, especially for A < Ac2, 
which may have led those authors to overestimate the 
gap in the paramagnetic state. 

Upon doping, two changes occur. First, the chemi- 
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TABLE I: Integrated density of states (both spins) for the 
(lowest- lying), "Zhang- Rice" (ZR) and "Upper Hubbard 
Band" (UHB) spectral features discussed in the text for U — 
9eV, A = 4eV (paramagnetic metal in undoped case) and 
A = 2eV (paramagnetic insulator in undoped case) at dopings 
indicated. 

cal potential moves into the Zhang-Rice band (hole dop- 
ing) or the upper Hubbard band (electron doping). For 
A ~ Ac2 the associated changes in chemical potential 
are substantial: the Fermi level, measured relative to the 
non-bonding oxygen peak, shifts by almost 2eV; for the 
A ^ Aci the Fermi level changes rather less. Second, 
as can be seen by inspection of the Figure O and from 
Table |T1 the relative strengths of the different spectral 
features evolve. The issue of the number of states cre- 
ated by doping has received some attention in the litera- 
ture as a signature of "Mottness" We find that the 
changes in the spectrum do not have a universal dop- 
ing or interaction-strength dependence; however in the 
"charge-transfer insulating" regime A < Ac2 , each doped 
hole adds roughly two states to the Zhang- Rice band and 
one to the upper Hubbard band, while doping with elec- 
trons does essentially the opposite. In the paramagnetic 
metal case the changes in electronic structure are larger. 



C. Self-energy and velocity renormalization 

To further probe the particle-hole asymmetry we show 
in the upper panel of Fig. [3] the self-energy calculated on 
the Matsubara axis for 0.15 electron and hole doping at 
A = 4eV ~ Aci. The near-perfect agreement between 
the results of ED and CT-QMC calculations serves as 
a test of the reliability of our results. The lower panel 
focuses on the low frequency behavior, presenting the 
doping dependence of —dTt/duj\ui^o estimated from the 
values of lm'S,{iujn) at the lowest two Matsubara points 
of the ED calculation for both A < Ac2 (paramagnetic 
insulator) and A ~ Ad (paramagnetic metal). At very 
low doping the estimate becomes unreliable because the 
(very small) Fermi liquid scale is not easily resolved so we 
do not present results. We see again that the self-energy 
is systematically larger for hole doping than for electron 
doping, with the difference being more pronounced for 
the paramagnetic insulator case. 

However, the self-energy is not necessarily the most 
relevant measure of correlation strength. In a multiband 
model such as the one studied here, the Fermi surface p = 
Pf is given by the solution of — ReS(w) — H(p) = 

at w = 0. The degree to which d states participate in the 
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FIG. 4: Doping dependence of velocity renormalization 
^•y^barc inverse of self-energy derivative — 1/(1 — 
dTi/du). Upper panel: A = 2eV; lower panel A = 4eV. 



FIG. 3: Upper panel: comparison of CT-QMC and ED self- 
energies over a wide frequency range calculated for U — 9eV, 
A = 4eV ~ Aci at dopings indicated. Lower panel: dop- 
ing dependence of —d'E/du>\ui^o estimated from lowest two 
Matsubara points obtained from U — 9eV ED calculations. 



Fermi surface states \ipipF)) depends on the self-energy, 
which of course changes across the charge-transfer gap 
and with doping. By expanding for uj near and p near 
PF and defining V'^^''' dU/dpand Z = 1 - dRe'E/duj, 
we find that the physical quasiparticle velocity v* is given 
by 

and the bare Fermi velocity by the same equation but 
with Z = 1. The bare velocities (defined here in physical 
units by multiplying the result above by the lattice con- 
stant 3.8A) have only about 5% doping dependence on 
either the electron on the hole doped sides, but change 
substantially as one goes from electron to hole doping: 
at A = 4eV we have v'°''"' = 4.3eV - A for the 0.15 hole 
doped and 3.5eV — A for the 0.15 electron doped calcu- 



lation while at A = 2eV we have v^''''' = 4.2eV - A for 
0.15 hole doping and 3cV — A for 0.15 electron doping. 

We see from Fig. |4]that for the strong coupling (A = 
2eV) case the particle-hole asymmetry in S is slightly 
overcompensated by a difference in d character of the 
ground state wave function, so that the particle-hole 
asymmetry in the ratio is rather smaller in 

magnitude and of opposite sign, compared to that in 
ImE(itij„), while in the more weakly correlated case the 
wave function changes induce an asymmetry in velocity 
renormalization which, while small, is larger than that in 
ImE. 



D. Optical conductivity 

We have calculated the optical conductivity and have 
verified the results via comparison of the integral of our 
calculated conductivity to the independently calculated 
"kinetic energy" and to the average of the renormalized 
Fermi velocity over the Fermi surface. We also compared 
results obtained directly on the real axis from the ED 
calculation to results obtained by analytic continuation of 
the Matsubara axis self-energy obtained from CT-QMC 
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calculations using the methods of Ref. [l^. 

The upper panel of Fig. [5] shows the conductivity in 
the near gap region for a doping of one hole per CU-O2 
unit in the antiferromagnetic phase at A = 2eV, 4eV and 
4.5eV. As noted above, for parameters A ~ Ac2 the gap 
in the antiferromagnetic phase is rather larger than the 
experimentally measured value ~ 1.75eV32. On the other 
hand, the gap value determined from the antiferromag- 
netic phase of the A ^ Ad calculation is in reasonable 
agreement with data. This comparison places the ma- 
terials clearly on the metallic side of the paramagnetic- 
metal/charge- transfer-insulator phase diagram, in agree- 
ment with previous analysis based on the one-band Hub- 
bard model»i? However, it is important to note that 
the calculated conductivities in this frequency range are 
about a factor of 2 smaller in magnitude relative to ex- 
perimental data (note that a normalization error means 
that the conductivity results of Ref. 18 are too large by 
a factor of two). Some of the difference may arise from 
transitions to bands not included in the present calcu- 
lation, but it is possible also that the Peierls phase ar- 
guments omit important interband matrix elements even 
within the space of states we consider. This is an impor- 
tant issue for further study. 

The lower panels of Fig. [5] display the doping depen- 
dence of the conductivity, calculated in the paramagnetic 
phase for A — 4eV. Electron or hole doping adds opti- 
cal absorption strength (associated with motion of doped 
holes) at frequencies < 2eV. For both values of A the 
integrated absorption strength (up to 2eV) is found to 
be comparable for electron and hole doping, despite the 
differences in self-energy displayed in Fig. [S] In this fre- 
quency range, the optical matrix element is proportional 
to the Fermi velocity which depends on the renormalized 
Cu-0 energy difference set by the parameter A — Rel]((jj). 
The shifts in ReE as the chemical potential is tuned 
from hole to electron doping leads to a bare Fermi ve- 
locity which is larger on the hole-doped side than on the 
electron-doped side as noted in the previous sub-section 
to a change in the d-character of the wave function; these 
effects lead to changes in the optical matrix elements 
which compensate to a considerable degree for the change 
in self-energy. One of us previously arguedii that the 
experimentally observed similarity in low frequency op- 
tical absorption between electron and hole-doped mate- 
rials implied that the correlation strength was about the 
same for two systems; we see that this argument must be 
treated with caution. 

A more striking change is that hole-doping but not 
electron doping activates a strong transition at about 
4eV between the non-bonding oxygen band and the near- 
Fermi-surface states. This feature is not observed exper- 
imentally; indeed a comparison of Figures 7 and 8 of 
Ref. shows that in the range O ^ 4eV the optical ab- 
sorption in electron doped compounds is slightly larger 
than in hole-doped compounds. We note that in a more 
realistic three-band model with oxygen-oxygen hopping 
included the spectral weight in this peak would be spread 
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FIG. 5: Upper panel: optical conductivity in near-gap region 
calculated for antiferromagnetic phase of three-band model 
from QMC calculation at T = O.leV, U = 9eV, carrier den- 
sity of one hole per Cu02 unit (undoped case) and A values 
indicated. Parameters: A = 4.5eV, Ed = — 9.1eV; A = 4eV, 
£d = — 8.8eV; A = 2eV, Sd ~ — 7.8eV. Lower panel: dop- 
ing dependence of paramagnetic conductivity for A = 4eV, 
U = 9eV from ED caIculation(red dash-dotted lines: un- 
doped; blue long dashed lines: ±0.10 doped; magenta short 
dashed lines: ±0.18 doped) along with antiferromagnetic con- 
ductivity in undoped case from QMC calculation (black solid 
lines, parameters are the same as in the upper panel). ED 
parameters: 0.18 hole doping, Sd ~ — 7.7eV; 0.10 hole doping, 
Ed = — 7.9eV; undoped, Sd ~ — 8.8eV; 0.10 electron doping, 
Ed = — 9.6eV; 0.18 electron doping, Sd = — 9.9eV. 



over a wider energy range. Further experimental and the- 
oretical investigation of this issue is important. 



IV. COMPARISON TO ONE BAND MODEL 

A basic issue in the physics of the cuprates is the re- 
duction of the three-band model to an effective one-band 
model.^''^''* We test the reduction by comparing our re- 
sults at A « Ac2 to results obtained by applying the 
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FIG. 6: Comparison between the three-band model and the 
one-band Hubbard model at dopings indicated. Only frequen- 
cies and energies relevant to one-band model are shown. The 
behavior of three-band model at other frequencies and ener- 
gies are similar to Fig. [2] and the lower panel of Fig. [S] Up- 
per panel: Spectral functions. Lower panel: Optical conduc- 
tivities. The one-band model is computed at nearest neigh- 
bor hopping t = 0.37eV, UcS = 12t = 4.44eV, T ^ O.lt = 
0.037eV using CT-QMC solver. Parameters: 0.10 hole dop- 
ing: Sd = -2.7t = -l.OOeV; undoped: ea = -Gt = -2.22eV; 
0.10 electron doping: Ed = —9.3t = — 3.44eV. The three-band 
model is computed at [/ = 9eV, A = 2.5eV, T = using ED 
solver. Parameters: 0.10 hole doping: Sd ~ — 6.8eV; undoped: 
Ed = — 8.4eV; 0.10 electron doping: Ed = ~9.5eV. 



QMC solver to a one-band square-lattice Hubbard model 
at Ucff = 12t « Uc2 and t = 0.37eV chosen to reproduce 
the sphtting between the centroids of the upper and lower 
Hubbard bands. Fig. [H] shows the comparison of spectral 
functions and optical conductivities. To obtain spectra 
for the Hubbard model we used the analytic continuation 
procedure described in Ref. [H, while the three-band re- 
sults were obtained using the ED solver at parameters 
U = 9eV, A = 2.5eV and T = 0. AU calculations are 
done in the paramagnetic phase. 



The upper panel of Fig. [6] compares the spectral func- 
tions. One sees immediately that the one-band and three- 
band models provide a reasonably consistent account of 
the spectra within a few eV of the Fermi energy. Some 
differences are evident. From the upper panels one sees 
that the one-band model has a slightly larger bandwidth 
than the three-band model (the difference is most evi- 
dent for the half- filled calculation). This difference may 
be interpreted as indicating that the three-band model 
at A « Ac2 is equivalent to a one-band model at a UcS 
a bit greater than Uc2- 

The lower panel of Fig [6] compares the conductivi- 
ties. In the half filled case the two models give a quite 
consistent account of the absorption above the Mott- 
Hubbard/charge-transfer gap edge. Note that because 
the lowest-lying gap excitation is not optically active, the 
difference in gap values noted in the previous paragraph 
is not easy to see in this panel. The one-band and three- 
band calculations are done with different methods but as 
can be seen, the conductivities are similar and we have 
verified that the spectral weight (integral of a up to say 
leV) are similar for the two models. However, for hole 
doping a new feature in the three-band model appears at 
H. ~ 4.5eV. This is associated with transitions from the 
non-bonding oxygen band, which is of course not present 
in the one-band model. 



V. CONCLUSION 

In this paper we have used the single site dynamical 
mean field method to solve the three-band model believed 
to be relevant to the copper-oxide materials. The method 
includes the physics of Zhang- Rice singlets (doped holes 
reside largely on oxygen sites and have spin opposite to 
the copper spin), reproduces the characteristic features 
of the spectrum, and reveals (especially for parameters 
in the "charge-transfer insulator" regime) a particle-hole 
asymmetry in the self-energy which however is largely 
canceled by the difference in d-content of the near-Fermi 
surface states, leading to very similar velocity renormal- 
izations between electron and hole-doped compounds. 
One important implication of this finding is that (at least 
for the correlation strengths we have studied) the fact 
that doped holes are "Zhang-Rice singlets" while doped 
electrons are just conventional doubly occupied sites does 
not imply a significant difference in the physics and sug- 
gesting that a reduction to an effective one-band model 
may be reasonable. The issue of reduction to a one-band 
model has been considered by many previous authors; 
however the discussion has been largely been couched in 
terms related to reduction to t-J-like models. We prefer 
to directly compare spectral functions and conductivity. 
We found that for A ~ Ac2 (paramagnetic insulator) a 
one-band model gave a quantitatively accurate descrip- 
tion of the physics at scales below about 4.5eV (above 
this scale the effects of non-bonding oxygen states which 
are not included at all in the one-band model become vis- 
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ible). Similar results (not shown) are found for A < Ad. 
It is likely that for parameters much deeper in the insulat- 
ing phase (which we have not investigated) the "Zhang- 
Rice" effects may be more important and the reduction 
to a one-band model may be more problematic. 

A basic question in cuprate physics is the effective cor- 
relation strength governing the physics of the low-energy 
particles responsible for superconductivity. We addressed 
this question by calculating the gap in the insulating 
phase. The gap is of a "charge-transfer" rather than 
"Mott-Hubbard" nature, as stressed by many previous 
workers. Our comparison to the one-band model shows 
that the charge-transfer gap can be used to extract an 
effective U which, when used in a one-band model, re- 
produces the low energy physics reasonably well. In con- 
trast to Ref. we find that antiferromagnetism has a 
pronounced effect on the magnitude of the gap in the in- 
sulating state. We find that there is no reasonable way 
to obtain a gap of the physical scale (between 1.5 and 
2eV) if the materials are assumed to be paramagnetic in- 
sulating side of the phase diagram; rather they must be 
taken to be more moderately correlated. This conclusion 
is in agreement with previous results. "'^^ 

Three features of the calculation suggest potentially in- 
teresting directions for future research. First, the quasi- 
particle band structure depends on the renormalized d 
level energy e*^ = Ed + ReS(w — 0). This quantity 
has (especially for stronger correlations) a noticeable de- 
pendence on doping and on which side of the charge- 
transfer gap the materials are on. It also changes dra- 
matically between antiferromagnetic and paramagnetic 
states. Further investigation of the physical consequences 
of these changes would be useful. This might produce 
additional insight into the fundamental electronic struc- 
ture dichotomy in many transition metal oxides between 
the dramatic evidence for strongly correlated behavior in 
high energy spectroscopies and the more nearly band-like 
behavior of the lower energy excitations. 

Concerning the conductivity, as is seen most clearly 



in Fig. [6l the calculated insulating state conductivity 
is only a few hundred inverse-ohm- inverse-centimeters, 
rather smaller than the experimentally measured values 
^ 800— 1000r2~^cm~^. Some part of the discrepancy may 
arise from other bands, not considered in our calculation, 
but it may also indicate a failure of the Peierls-phase ap- 
proximation to the conductivity. A further issue in the 
comparison of the high-frequency calculated conductivity 
to data is the appearance, for hole doping but not elec- 
tron doping, of a strong feature relating to transitions 
from the non-bonding oxygen bands to the near-Fermi- 
surface states. This feature implies that the high energy 
conductivity in the 4-6eV range should be greater for hole 
doped than for electron doped compounds. This is not 
seen experimentally. Inclusion of oxygen-oxygen hopping 
will spread the excess spectral weight over a wider fre- 
quency range, perhaps mitigating the discrepancy with 
experiment. Further, one must bear in mind that the 
electron and hole doped materials studied experimentally 
have different crystal structures, introducing further un- 
certainties in the comparison. However, the issue war- 
rants further study. 
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